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We derive a low-energy quantum field theory from quantum chromodynamics (QCD) that holds 
in the limit of a very large coupling. All the parameters of the bare theory are fixed through 
QCD. Low-energy limit is obtained through a mapping theorem between massless quartic scalar 
field theory and Yang-Mills theory. One gets a Yukawa theory that, in the same limit of strong 
coupling, reduces to a Nambu-Jona-Lasinio model with a current-current coupling with scalar-like 
excitations arising from Yang-Mills degrees of freedom. A current-current expansion in the strong 
coupling limit yields a fully integrated generating functional that, neglecting quark-quark current 
coupling, describes all processes involving glue excitations and quark. Some processes are analyzed 
and we are able to show consistency of Narison-Veneziano sum rules. Width of the a resonance is 
computed. The decay r\ — > rj + tv + + tv~ is discussed in this approximation and analyzed through 
the more elementary processes rj — >• r\ + a and a — > ir + + 7v~ . In this way we get an estimation of 
i-O ■ the mass of the a resonance and the value of the rj decay constant. This rj' decay appears a possible 

' source of study for the o resonance. 
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I. INTRODUCTION 

^ ; 

'' With the discovery of asymptotic freedom 0, Q and the settlement of a theory of strong interactions with quantum 
CV chromodynamics (QCD), the question of the determination of the low-energy behavior of the theory has been increas- 
ingly demanding. Once it was clear that the high-energy behavior of strong interactions was properly described, the 
emergence of the low-energy spectrum of the theory must have been obtained. This question encountered a serious 
difficulty as perturbation theory is useless in this case due to the strength of the coupling and there were no useful 
mathematical techniques at that time. At the very start of this story, Wilson come out with a discretized approach 
for QCD Q . This approach made the basis for the birth of lattice QCD that still today is the most reliable approach 
to solve QCD at low energies. 

In the course of time, some approaches have been devised to manage low-energy QCD. The idea behind this is 
to move to a phenomenological theory having all the properties of QCD. With this view in mind, a very successful 
approach was devised by Weinberg Gasser and Leutwyler Chiral perturbation theory (ChPT). Since then, this 
approach has found increasing success (for a review see |y|) and experimental support Q- This approach introduces 
some phenomenological constants that are fixed through experimental measurements. Analyses of weak hadronic 
T — I , decays are accomplished through this approach [8] . 

On a different side, Nambu-Jona-Lasinio model [{J [l(| proved to be a very satisfactory description of hadronic 
phenomenology [TTI j . In this case the model is based on a quark current-current interaction and is phenomenologically 
built to display dynamical chiral symmetry breaking. In its original formulation, Nambu-Jona-Lasinio model is neither 
renormalizable nor confining. A proposal to modify this situation has been given in [f2| and further generalized with 
a proper gluon propagator and starting from QCD in [1 31) . On a Nambu-Jona-Lasinio model bosonizationc can be 
applied to introduce meson fields from quark fields [Tj]. Bosonization techniques make comparable Nambu-Jona- 
Lasinio models and ChPT [l5j . 

Not less relevant has been the approach of quantum spectral sum rules that have taken the start from the pioneer 
work of Vanshtein, Shifman and Zakharov [l6|,[l3] that also pointed out the existence of gluon and quark condensates. 
Sum rule techniques showed the need for a further meson with a mass around a GeV [l8j. Bosonization procedures 
produce the appearance of the field of the a meson in the Nambu-Jona-Lasinio model and this field is needed in ChPT. 
Indeed, existence of this particle has been the argument of longstanding discussion and, after an initial insertion, it 
was dismissed by Particle Data Group for over twenty years. Problem in experiments arisen from the fact that this 
resonance is quite broad and so difficult to be properly measured. Today there is widely acceptance and a decisive 
step forward is due to the work of Caprini, Colangelo and Leutwyler (i~9| further supported by the work of Yndurain, 
Garcia-Martin and Pelaez [2(ij |. The main question now moved from the existence of this resonance to its nature. In 
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this paper we show evidence toward a large gluonic component in it. 

Notwithstanding a wide success encountered by these techniques, the understanding of the lower part of the spectrum 
of QCD escaped us and is reason of a very hot debate in these days [21,]. The reason has been already pointed above 
and relies on our difficulties to manage QCD at very low energies. Lattice computations have given very few help 
about and it is possible that there is something deep to be better understood with this approach. 

So far, a fundamental question remains unanswered. There has been no approach to derive a low-energy limit 
directly from QCD. There has been no successful and widely accepted proposal so far. But a significant track has 
been unveiled by Goldman and Haymaker (22|. These authors proved that QCD has a meaningful low-energy limit 
provided the form of the gluon propagator is given and is not as singular as the propagator of a free massless particle 
in the limit of lower momenta. On a same track, one can do a current expansion into the generating functional of QCD 
recovering Goldman and Haymaker setback [23|. At some point one should know the gluon propagator to recover a 
low-energy limit from QCD. 

The question of the understanding of the gluon propagator in the low-energy limit has been matter of wide debate 
in the last few years, for a pure Yang-Mills theory. There has been groundbreaking theoretical work, well recounted in 
[24j . claiming that, lowering momenta, the gluon propagator should have gone to zero. Similarly, as a consequence, the 
running coupling should have reached a fixed point. People working on the lattice, at first, were not able to have a clear 
understanding about the right behavior and, for some years, appeared as this theoretical scenario was supported. A 
crucial step happened at the XXV International Symposium on Lattice Field Theory, held at Regensburg in Germany 
on 2007, where two groups announced results on very large lattices |25|,|26|. Results on lattice showed unequivocally 
that the gluon propagator, at very low momenta, reaches a finite value different from zero, with a curve resting on 
a plateau for the Landau gauge. Running coupling is seen to tend toward zero having no fixed point at all. This 
latter result is in agreement with the beta function for a SUSY Yang-Mills theory that is exactly known [13] and its 
extrapolation to Yang-Mills theory [2^, [2!| . 

In this paper we will derive the low-energy limit of QCD, based on a sound derivation of the gluon propagator in 
the infrared. Several applications to low-energy QCD processes will be given. The paper is so structured. In seclTIl 
we introduce our notations for QCD. In sec lllll we perform the analysis of the low-energy limit of QCD obtaining the 
gluon propagator and showing how to overcome Goldman and Haymaker setback. In sec II VI we gives some applications 
to low-energy processes in QCD where glue excitations and quarks interact at tree level. Finally, in secfVl conclusions 
are presented. 



II. QUANTUM CHROMODYNAMICS 



In order to fix notation, we give below the formulation of QCD we will work with. 
The action for QCD can be written as 



~>YM 



d 4 x 



being 



^G^ + l-id-A) 2 



(1) 
(2) 



— — ®i>A* + gf a "ApAy 

the field tensor and g the coupling constant that is dimensionless in four dimensions. f abc are the structure constants 
of SU(3) group with the generators [A a ,A{,] = if abc \ c . Here a,b,c— 1,2,3. Finally, the value of a gives the fixing 
gauge term. 

For quarks one has 



S m = J d^x Q (ij ■ 9 



+ g^-l-A a -m q ) q. 



being q = u, d, s, . . . the quark fields and m q the corresponding masses. 

Finally, we have to introduce Fadeev-Popov ghost. This is done in the following way 

S g = - J d*x [d^d^d 1 + gf abc d^A htl c c ] . 

From these equations we can formulate the generating functional for quantum field theory as 



Z[rt, fj, e, e, j] = J JJ [dq] [dq] [dc] [dc] [dA] , 



,i(S YM +S m +S g ), 



(3) 



(4) 



(5) 



Q 



e ifd xJ2 q k(x)n q (x)+V q (x)q(x)] e iJd 4 -xlc a (x)e a (x)+e a (x)c a (x)] e ifd' l xj? l (x)A a > J -(x)^ 
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here and in the following h = c = 1. Our aim in this paper will be to get a proper characterization of the behavior of 
the quantum theory in the low-energy limit. 



III. ANALYSIS OF LOW-ENERGY LIMIT 



A. General considerations 



Quantum field theory relies heavily on our ability to find out exact solutions to classical equations of motion. 
Once such solutions are known, we are able to do perturbation theory and obtaining results to be compared with 
experiments. It is interesting to note that the kind of perturbation series will depend on the nature of the exact 
solutions we will work with. It is also possible to work with approximate solutions as starting point for perturbation 
theory. In this case one must be sure to work consistently and be able to get systematically all higher order terms. 

A possibility to formulate a strong coupling expansion in quantum field theory was put forward in [30]. In this 
case we were able to apply this approach to a quartic scalar field theory obtaining the corresponding pr opagat or and 
the mass spectrum. The theory is seen to be trivial as, at lower momenta, the coupling goes to zero [3ll433| . The 
corresponding expansion is indeed a strong coupling expansion as the parameter of the series is the inverse of the 
coupling. 

Due to the peculiar nature of a quartic field theory and a simple observation that Yang-Mills theory seems a 
vectorial generalization of this theory, one may think that a relation should exist between these twos. Indeed, there 
exists a class of exact solutions that are common to the classical equations of these theories but these solutions can 
only depend on time variable. When a gradient expansion is used to correct these solutions one observe that this 
mapping is only perturbative in the inverse of the couplings. Couplings are connected by the relation A = Ng 2 being 
A the coupling of the scalar theory and g the coupling of the SU(N) Yang-Mills theory. This is the content of the 
so-called mapping theorem that we were able to prove in [34l . |35| . The relevance of this result relies on the fact that we 
can pass all the results of the scalar field theory to the Yang-Mills theory when a strong coupling limit is considered. 
In this way we are able to get both the gluon propagator and the mass spectrum in the same limit and Yang-Mills 
theory is seen to display a mass gap. 

An immediate by-product of the mapping theorem is the opportunity, in the vein of Goldman and Haymaker [22| . 
to obtain the correct low-energy limit of QCD with all the parameters properly fixed by the theory. This limit is a 
Nambu-Jona-Lasinio model equivalent to a Yukawa model in a proper approximation. In the next sections we will 
exploit all this in a consistent mathematical way. Our main goal will be to put on a sound basis the techniques to 
analyze in a closed form low-energy QCD. 



B. Exact solutions and mapping theorem 

Quite recently, we showed that a gradient expansion is indeed a strong coupling expansion [36l|. This can be seen 
very easily for a quartic scalar field theory. One has 

□</> + A0 3 = (6) 
with A — > oo. Let us rescale the time variable as r = s/\t and take 

' ofiV (7) 



A I A 2 



dfyo + 4>l = (8) 



We get a consistent set of differential equations 

9^ + 3^1 = A 2 o (9) 

(10) 

that is clearly a gradient expansion. This means that gradient expansion is nothing else than a strong coupling 
expansion [36]. Being this a gradient expansion, it is not difficult to note that an exact solution to our equation of 
motion is obtained if only time dependence is retained solving just the first equation of the given set. Moreover, this 
will produce an exact solution to the full equation after a Lorentz boost. Such a solution can be written down as 

<P(t,0)=[iCj ) sn ( ( - ) /»/ hH.i J 
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being /i and 9 integration constants. A Lorentz boost takes this solution to 

<f>(x)= f i(^) Z sn(p-x + 6,i) (12) 

i 

provided p 2 = (§) 2 p 2 ■ Now, given Yang-Mills equations 

d^d^,A a v - (l - &,(0M£) + g] ahc A^id^Al - d v A c ^ + gf ahc d^{A h ll Al) + g 2 f abc f cde A b > i A d l Al = (13) 

it is always possible to choice a set of constants 77° , Smilga's choice [37} , such that the solution (fTTj) is also a solution 
of these equations. But here we cannot just boost these solutions to get something like (fT2"j) . Rather, due to the 
presence of an arbitrary choice of gauge, we can only state the following perturbative mapping between the scalar 
field and Yang-Mills potentials 

AZ(x)=r,Z<t>(x) + o( n L-) (14) 



Ng 

provided A = Ng 2 . This is exactly the content of the mapping theorem [34|, [35[. Due to the nature of these 
approximate solutions of the Yang-Mills equations, they indeed provide the behavior of the classical theory when the 
coupling becomes increasingly large. This is exactly the limit we are interested in. So, a quantum field theory can be 
built on these solutions through a perturbation theory in the strong coupling limit. 

C. Gluon propagator 

In order to work out a low-energy limit for QCD we will apply the mapping theorem to the generating functional 
(0). This will give 

ZhfjJ] = I nW^W]e pM,/A ^^^^]x (15) 



e if r^zE, q(x)[l\id+g^4:)-m q ]q(x) e i j d i xY. q k( x )Vq(x)+ , nq(x)q(x)} x 

having set = r\ ■ j . It is interesting to note that, at this order of the strong coupling expansion, the ghost field 
decouples and can be ignored in our analysis, this is in agreement with lattice results j25l |26|. Higher order terms can 
be properly computed through a gradient expansion of the classical equations of motion. We can recognize here the 
quantum field theory of a Yukawa model. This model can be cast into a Gaussian form and so can be reduced to a 
manageable form. So, let us consider the generating functional for the scalar field 

Z \j] = J ^Je^MlW) 2 -^^] (16) 

and take the limit A — > 00. For a strong coupling expansion, we separate gradients from other terms as 

Z \j] - / dW'^H^f-^+^e-^^i^) 2 ] (17) 



and finally, we have to solve 

d 2 <j> c + Ml = j. (is) 

Taking the expansion 

= cb c + Scj) (19) 

the generating functional is reduced to the form [3(| 

Zo[j] =e^ d ^ F[j] (20) 
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with the functional F[j] containing terms negligible in the limit A — > oo. Another approximation to be applied is to 
take, in a small time limit 38, 391, 



(21) 

(22) 
(23) 

(24) 



M x )~ J d 4 x'G(x-x')j(x'), 

also to be seen as a current expansion [23^. being 

G(x-x') = S 3 (x-x')G{t-t') 

and 

d^Git - t') + XG(t - t' f = 6{t - t') 

that can be exactly solved giving 

G(t - = 6{t - t> (f ) 4 ™ C 4 a**, * 



and the corresponding time-reversed solution. This gives finally a Gaussian generating functional that holds for a 
strongly coupled theory, the infrared limit, 



and the propagator given by 



Z [j] w S di ' xdi y3{x)&(x~y)i{y) 



A(a; - y) = <5 3 (a; - y) G+(t x - ^) + G-{t x - t y ) 



(25) 



(26) 



with G+ given in eq.tfMf and G_ the corresponding time reversed, t — > — t, Green function. Now, we are able to 
obtain the complete propagator of this quantum field theory in the infrared limit. The relevance of this result relies 
on the fact that, due to the mapping theorem, this will also represent the infrared gluon propagator for Yang-Mills 
theory. Indeed, we note that for the Jacobi snoidal elliptical function the following identity holds [40] 



, . 27T ^ (-l)n e -{n+^ 
SmuA) — — > -77 — sin 

v ' K(i) ^ 1 + e -( 2n + 1 ^ 



(2n+ 1) 



2K(i) 



being K (i) the constant 



So, one has 



K(i) 



d6 



o V1 + sin 2 



A( W ,0) = £ 



1.3111028777. 



n=0 



ui 2 — mi + ie 



being 



and, 



Br, 



7T 2 (-l)«+l e -(»+!K 

(2n + l >K 2 {i) l + e -(2n+i> ' 



m n — \ n 



1 \ 7T /A 



2 / JT(i) \2 



(27) 



(28) 



(29) 



(30) 



(31) 



the spectrum of the theory. This theory displays a mass gap in the infrared limit and this also must happen to Yang- 
Mills theory. Now, we can boost this propagator from the rest frame and this will amount to sum all the gradient 
terms providing in the end 



A(p) = £ 



n=0 



B n 

p 2 — m 2 x + ie 



(32) 



6 



So, we are in a position to write down the infrared propagator for SU(N) Yang-Mills theory in the Landau gauge as 

Df v {x -y) = 5 ab (r,^ - Mnj A(p) (33) 

provided we take A = Ng 2 . This is a key result, as it overcomes Goldman and Haymaker setback, and we are able to 
get a meaningful low energy limit for QCD. 

D. Low-energy QCD 

Once we know the gluon propagator, we can do a change of coordinates into the generating functional of QCD. 
This change of coordinates is the one suggested in the analysis of the scalar theory. So, in the Yukawa model (p~5)) we 
put 



x) ps / d 4 yA(x - y) 



Q 



(34) 



and we use the reduction to a Gaussian form exploited in the preceding section. This will give the generating functional 

Z[r),fj,j] = e i{N2 - 1) f dixdiv ii j ^ x)A(x - y)j *W} f Y[[dq][dq]x (35) 

£ i J ■d*xj: 9 q(x)[y.(id+g?%L J 'AA(i-»)[i*(9)+|9E,?l!/)^5(ri])-m t ]?(»:) e i/'i 4 *E,[?W') J W+S»W?(*)] x 

+0(l/VNg). 

This result is quite striking. The Yukawa model we obtained through the mapping theorem just gives rise to a 
modified Nambu-Jona-Lasinio model, without contact interaction, plus a term describing the pure interaction of 
the quark fields with a set of scalar excitations. These excitations should represent the possible observed glueball 
spectrum. For completeness, we give below the Lagrangian of the model, neglecting the Gaussian term describing 
free scalar excitations, 

L = ^2q(x)(ij- d-rn q )q(x) (36) 
i 
9 



2 



^q(x)X a V a ^q(x) J d 4 yA(x - y) H {y) 



+ JtE*) A "</9W / d 4 yA(x- y )q>(y)\ b r,lYq'(y) 



2 4 

9,9' 



+ 0(l/VNg) 

This Lagrangian, together with the given generating functional, fixes univocally the behavior of QCD at lower energies. 
It is interesting to note that we have no yet done a Nambu-Jona-Lasinio point contact approximation. So, we can be 
sure that the theory is renormalizable. The form of the gluon propagator grants that a point contact approximation can 
always be done recovering a fully Nambu-Jona-Lasinio model removing Goldman and Haymaker setback. Confinement 
is granted again by the gluon propagator that displays a full set of massive excitations. When we do a reduction to 
the Nambu-Jona-Lasinio model with a point contact interaction we lose both renormalizability and confinement but 
breaking of the chiral symmetry is maintained. 



E. Nambu-Jona-Lasinio limit 



Nambu-Jona-Lasinio limit corresponds to current-current approximation in the low-energy field theory we obtained 
so far. Existence of this limit can be easily inferred from the gluon propagator (|32[) that has a finite limit when 
momenta go to zero. One has 

3 7fi 

A(x - y) « — 8\x - y) (37) 
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being a the string tension as defined above. From lattice computations one has ^/a — 0.41 -j- 0.44 GeV but this value 
should be obtained from experiment. Now, we define 



G 



NJL 



and 

The corresponding Lagrangian will be 

Lnjl = 



3.76 



3.76 



3.76- 



4:77 a „ 



g 'njl = - — 9 = 3 - 76 



. V47ral 



^2q(x) (i'y ■d-m q )q(x) 
i 

lG' NJL J2Q(x)* a v a ^q(x)U(x) 

q 

IGj^l ^ q {x)X a v a^ q{x) q' {x)X ^q'(x) 



(38) 



(39) 



(40) 



q.q' 

+ 0(l/VNg) 



producing a generating functional 



Z NJL [r),fj,j] = e <^-VJ** d HhM*)M*-v)Mv) 



Y[[dq}[dq] 



(41) 



e if d*x [E« g(*)(i7-a-m g )g(x)+AG' NJ£ E, ?(*)A <, <T / *g^)jV(x)+i^ i E g , 9 ' «(*) A^VM*)?' (x) A\> 7 V (x)] x 

At this stage we take a current expansion [23j as igj^, J, and so on and we are left with a first order approximation 
in the strong coupling limit 

Z [r,,fj,j] « e^-D/^liwWA^M] [l[[dq}[dq}x 



(42) 



e i / d 4 a: E 5 [g(x)j7 9 (x)+^(x)?(a;)] _ 

where only the term jgj^, is retained. This approximation amounts to ignore quark loops and to consider only 
quark-gluon vertexes. We can perform the integral if we know the solution to the equation 



1 



ll .d-m q + -G' NJL X a V a ^ U (x) 



S q [j<j>;x - x] = 8 A (x - x'). 



(43) 



The idea at this stage is to perform again a gradient expansion retaining just the leading term. One gets without 
difficulty the following result 



A" 



iG'NjL—lQlir)'} I dtfj^t' -t x ,,x-x') 



(44) 



S\j<t>,x- x'] = 9(t x - t x ,)5 3 (x - x')exp 
Finally, we have the low-energy generating functional for QCD as 

Zo[v,V,3</>] ~ exp |^(^ 2 - 1) J d 4 xd 4 yj 4> {x)A(x - y)j<f,(y)^ exp ji J d 4 xd 4 y^ Vq(x)S[j 4> , x - . (45) 

This generating functional describes the behavior of glue scalar excitations and their interaction with quarks. To 
account for loop contributions due to quarks and scalar excitations one has to consider again quark current-current 
interaction. 
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IV. QCD PROCESSES AT LOWER ENERGIES 
A. Narison-Veneziano sum rules 

In the eighties, Narison and Veneziano derived some sum rules for scalar resonances in QCD [4l|. These rules are 
relevant to obtain some properties, as the decay widths, of such particles and give hints for their observation. In the 
following we prove that Narison-Veneziano sum rules are properly verified for the spectrum of ++ mesons decaying 
in 7r + 7r~ obtaining in this way some relevant information about this decay but, more importantly, we get the form 
factors of scalar resonances. 

In order to obtain this result we use the leading order approximation of an expansion in \j\[Zg for QCD. Firstly, 
we consider 

jE^+--^r = 1 (46) 

s s 

and from this we get the corresponding value of a s . We use this to finally check that 

J V2/ 5 «0. (47) 

s 

In order to prove this, we have to compute fg that are the decay form factors for scalar mesons. This is easily 
accomplished through the mapping theorem that gives 

A a M (x)=ffc<Kx) + 0(l/V3g) (48) 

being 

□0 + 3 ff 2 = O. (49) 

We recognize here the one-point function of the quantum field theory when all the currents are set to zero. This 
equation has the following exact solution 



cj)(x) = J-^-sn{p-x + e,i) (50) 
being sn a Jacobi elliptic function and provided that 

p 2 = a (51) 

being \fd = A(3<? 2 /2)3 the string tension taken to be 410 to 440 MeV in lattice computations. But this solution can 
be also given as the sum of free massive excitations through its Fourier series as 

I ^ 9-7T ,= ~(™+3) 7r 

*'°> = fern B-i)"^^— • + * + < 52 > 



n=0 



with the meson spectrum (see (|31|l) 



and the form factors 



M s „ = (2n+l)^yV^ (53) 



/s " = V 67ra s tf(i) ( " 1) "l + e-( 2n + 1 >' r (M) 

where we see that they depend on 1/ ^/a^ that is clearly not obtainable from a standard small perturbation approach. 
It is interesting to note that they have alternating signs in agreement with what one should expect from Narison- 
Veneziano sum rules. 
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Finally, putting all together and doing the severe approximation that all gs-x+n- are equal, evaluating the sum will 
give from eq. pS)) 

9S " + *-* (6 ™ s)h ^ki)^ (55) 
that, taking gs-n+n- ~ 3.8 -j- 3.9GeV, will give 

a s = 0.25 ^0.26 (56) 
a fairly good value. To complete our check, we consider the other Narison-Veneziano sum rule to obtain 

9s*+k- .^.x 0.024 « 0.02 GeV 2 (57) 

(67TQ! S )2 

and so we see that both sum rules are well verified for ++ glue excitations. 

B. Properties of the a meson and higher scalar excitations 

Several processes involving glue scalar excitations can be computed through the generating functional (|4"S")) . The 
mapping theorem showed that the glue sector of the QCD action produces a set of scalar excitations with a mass 
spectrum 

mn = [2n + l) W{i)^ (58) 

being 

a = V67ra s /x 2 (59) 

the string tension. From these considerations, we observe that the spectrum displays a massive state with a mass 
around 500 MeV. Such a state has been observed and some phenomenological analysis strongly support its existence 
[l9l I20I ] . This state is identified as the a meson. This resonance is broad and has a prevailing decaying mode in two 
pions and a very small decay rate in two photons. This latter decay has been properly understood as a rescattering 



effect 42]. In the debate about the nature of this resonance there have been some works pointing out its mostly glue 
content |42l - l44| . We will give a significant support to these results. 
The process 

a -> 7T+ + n~ (60) 
can be described by the Feynman diagrams given in fig HI This is a vertex process and the amplitudes are given by 

S 5 id 6 iS 



5j$(x) i$r) u (x) 8r)d(x) iSr)d(x) 5rj u (x) 
X a 

i-G'NjL—lwnt x 



(61) 



[9(t 4 - t 3 )6(t 2 - t 5 )9(t 2 - t 6 - t 1 )5 3 (x 2 - x 5 )6 3 (x 4 - x 3 )S 3 ( Xl ) 
+ 6{U - h)9{t 2 - t 5 )0(t 4 -t 3 - h)S 3 (x 2 - x 5 )6 3 {x 4 - x 3 )S 3 ( Xl )] . 

In order to evaluate this vertex we use the following expression for the Heaviside function 

1 r + °° 1 

Bit) = / dE e- lEt . (62) 

We have to go from this n-point function to the probability amplitude introducing pion fields. We do this with the 
following rule: we remove the contributions from the Heaviside functions with the product m a f%, assuming a at rest 
and m a its mass and the pion decay constant taken to be 130 MeV. This corresponds to LSZ reduction. So, when 
we take the trace of the square of the amplitude, we will be left with 

\M if \ 2 = A{N 2 -l)G' 2 NJL mlJt (63) 
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FIG. 1: a — > 7r + + tv amplitudes. 



that we specialize to N = 3. This will give the rate 



_ 2 ^/2 

— — ^ JVJL 
7T 



Ami 



"IT, 



(64) 



that is in agreement with recent derivations from experiments [2(| being T a /2 = 255 ± 10 MeV for a mass m CT = 
484 ± 17 JVfeV when a s 0.59. A similar conclusion can be drawn also with respect to the derivation given in 
(l9j . presenting a mass vn a — 441^g 6 MeV and T a /2 = 279j]^ 2 5 MeV, for a s w 0.47 showing that both results arc 
consistent each other with respect to physical values of the strong coupling constant. In these computations we have 
fixed the string tension through the mass formula for the a resonance obtained putting n = into eq.(|3ip. 
From this result, we are able to evaluate the coupling as given in literature [42|, S3]. One sets [HI 



r„ = 



16nm 



4 m 2 



(65) 



being g^ir the coupling. This gives 



AV2G' NJL m a fl « 1.68 4-1.93 GeV 



(66) 



when the parameter given in [19j, [20J are used. As one should expect, this number is quite sensible to the value of the 
a mass and the running coupling. 
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In order to have a comparison with the observed spectrum of Yang-Mills theory on the lattice [45|, |46| , we consider 
the standard two-point function 

Cf l/ (t-t') = (A-(t, 0)4(^,0)) (67) 

having used translational invariance. Now, we apply the mapping theorem and we have in the end 

C;l(t-t') = riXm,0mt',0)) + O(l/VN9). (68) 

This result is rather beautiful as we can use the exact classical solutions of the scalar field. Indeed, this function 
solves the equation df((f>(t, 0)<fi(t' , 0)) -I- \(4>(t, 0)(f>(t' , 0)) 3 = S(t — t') and so we can conclude that 

C;l(t - t>) = 6(t t' Kv b v E(- 1 )" ^ ( , )l + e -( 2 n +1 ), e ' tm " (t ~ + c - + 0(l/VNg) (69) 
as expected from lattice computations where the spectrum is given by 

1 

being yfa = fJ- the root of string tension. So, the mapping theorem grants an identical excitation spectrum 

for a pure Yang-Mills theory and a quartic massless scalar field theory in the infrared. This has been confirmed quite 
recently on lattice computations for d=2+l 47|. In d=3+l, for comparison with (45[, we get the following table 



Excitation 


Lattice 


Theoretical 


Error 


a 




1.198140235 




0++ 


3.55(7) 


3.594420705 


1% 


0++* 


5.69(10) 


5.990701175 


5% 


a* 




2.396280470 




2++ 


4.78(9) 


4.792560940 


0.2% 


2++* 




7.188841410 





where our numbers are computed with the formulas 



m n ( 1 



V 2; K(i) 
[n + m + 1) 



(71) 



fa v y K(i) 



that arc pure numbers (see also (48[). 

Finally, we note that a decay S„ — > tt + + 7r _ is a clear signature for these scalar excitations. This should explain 
why it is so difficult to identify them clearly with respect to other typical processes in strong interactions. It should 
also be noted that mixing may be at work worsening the situation. 



C. 7/ decay 

Computation of the decay of 7/ to r]ir + ir~ is an interesting question in QCD. Chiral perturbation theory, at the 
leading order, is very far off the measured value [4{|. So, it is interesting to see how low-energy limit of QCD recovers 
the amplitude. In our case, the diagram is given in fig|2]In this case there is a virtual a meson exchanged and this is 
a two-vertex process. In the following, using our approach, we exploit the possibility that 7/ is emitting a a meson 
and we get the rate. Indeed, we can cut this diagram in two simpler ones in figEl From this latter diagram we 
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can consider rf decay just as a production of a particles and we can immediately put down the amplitude for the 
production process rf — > rj + a as 

|A%| 2 = (iV 2 - \)G'\ JL m%,f a f n (72) 
remembering that now this is a two-body decay. So, the width will be given by 



1 2 / m\, + m4 + ml - 2m 2 ,77j2 - 2m 2 ,, m 2 , - 2m\ml 
n ' ~ 2tt NJL ' m " n '1° h Y ^ ' ^ 



From eg. (1541) we have 



cr 27T e 2 

^-^^fe^yiT^ (74) 



and we are in a position to give an estimation of the mass of the a resonance assuming that the threshold must be 
overcome for this process to be seen as is. Using the mass formula for the a we get the limit values yja — 0.342 GeV 
being m„ = 0.410 GeV that are quite near current estimations (e.g. see [13]). Finally, taking f v ps 0.019 GeV, very 
near f„i w 0.024 GeU [50| . and m„/ = 0.958 GeV we have in the end T ss 91 fceU in agreement with the experimental 
value (see [5l|). This result is quite satisfactory and could give a proper account of the observed process. Indeed, with 
our understanding of this decay, we will have r^/ W7r = T v iT a /(T v i + r CT ) w r,,/ taking e.g. Tcr/2 = 279l^ 2 5 MeU as 
given in [ll|. We can conclude that this process could be quite effective to extract the properties of the a resonance. 
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V. CONCLUSIONS 



We have given a consistent low-energy limit of QCD when a strong coupling approximation is considered. In this 
paper we have shown simple applications to work out the consequences of our approximations. This technique could 
be extended at higher orders but having a well-definite model at the leading order can already help to work out most 
of the phenomenology observed in current experiments. 
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